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Abstract. The NA57 experiment has collected high statistics, high purity samples 
of Kg and A, 5 and SI hyperons produced in Pb-Pb collisions at 158 A GeV/c. In this 
paper we present a study of the transverse mass spectra of these particles for a sample of 
events corresponding to the most central 53% of the inelastic Pb-Pb cross-section. We 
analyse the transverse mass distributions in the framework of the blast-wave model for 
the full sample and, for the first time at the SPS, as a function of the event centrality. 
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1. Introduction 

The study of ultrarelativistic heavy-ion collisions is motivated mainly by the quantum 
chromodynamics (QCD) prediction that at sufficiently high energy density the excited 
nuclear matter undergoes a phase transition into a system of deconfined quarks and 
gluons (quark-gluon plasma, QGP) PQ. 

Strange particles have proved over the past years to be a powerful tool for the study 
of reaction dynamics in high-energy heavy-ion collisions. 

The WA97 experiment has measured |2] an enhanced production of particles carrying 
one, two or three units of strangeness in central Pb-Pb collisions at 158 A GeV/c with 
respect to proton induced reactions (strangeness enhancement). The observed pattern 
of the enhancements as a function of strangeness was predicted more than 20 years ago 
to be a consequence of a QGP formation [3]. 

NA57 is a dedicated experiment at the CERN SPS for the study of the production 
of strange and multi-strange particles in Pb-Pb collisions [4 . It continues and extends 
the study initiated by its predecessor WA97, by (i) enlarging the triggered fraction of the 
inelastic cross-section thus extending the centrality range towards less central collisions 
and (ii) collecting data also at lower (40 A GeV/c) beam momentum in order to study 
the energy dependence of the strangeness enhancements. 

In this paper we concentrate on the analysis of the transverse mass (m T = 
a/]?t + m 2 ) spectra for A, E~, Vt~ hyperons, their antiparticles and Kg measured in 
Pb-Pb collisions at 158 A GeV/c. 

The m T spectra are expected to be sensitive to the details of the production 
dynamics El- I n particular, for a fireball in local thermal equilibrium, the shapes 
of the m T spectra depend both on the thermal motion of the particles and on the 
collective flow driven by the pressure. To disentangle the two contributions, namely 
thermal motion and transverse flow, one has to rely on models. Several freeze-out 
scenarios have been proposed (for an overview refer to 7j). They differ in the freeze- 
out geometry and in the flow velocity profile. In the present analysis we consider the 
blast-wave model which assumes cylindrical symmetry for an expanding fireball 

in local thermal equilibrium, with different hypotheses on the transverse flow profile. 

The paper is organized as follows. The NA57 apparatus is briefly described in 
section 2. Section 3 deals with the measurements of the transverse mass spectra. 
Inverse slope parameters are presented and discussed in section 4. The blast-wave model 
is introduced in section 5 and compared to data under different conditions. Finally, 
conclusions are drawn in section 6. 

2. The NA57 experiment 

The NA57 apparatus, schematically shown in figure HI has been described in detail 
elsewhere jHJ- The A, S~, Q~ hyperons, their antiparticles and the Kg mesons are 
identified by reconstructing their weak decays into final states containing only charged 
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Figure 1. The NA57 apparatus placed inside the 1.4 T field of the GOLIATH magnet. 

particles: 

K s 7r+ + n~ A -* V + ^~ 

5" -> A + 7T- Q- ->• A + K- ^ 

u 

and the corresponding charge conjugates for antihyperons. 

Tracks are measured in the silicon tracking telescope, an array of pixel detector planes 
with 5x5 cm 2 cross-section, with a total length of about 30 cm. To improve the 
momentum resolution of high momentum tracks a lever arm detector -an array 
of additional pixel planes and double-sided silicon microstrip detectors - - is placed 
downstream of the tracking telescope. The whole silicon telescope is placed in a 1.4 
Tesla magnetic field, above the beam line, inclined and aligned with the lower edge of 
the detectors laying on a line pointing back to the target. The inclination angle a with 
respect to the beam line and the distance d of the first pixel plane from the target are 
set to accept particles produced in about half a unit of rapidity around central rapidity 
and medium transverse momentum: a = 40 mrad and d — 60 cm. 
An array of scintillation counters (Petals), placed 10 cm downstream of the target, 
provides a fast signal to trigger on the centrality of the collisions. The triggered fraction 
of the Pb-Pb inelastic cross-section is about 60%. The centrality of the Pb-Pb collisions 
is determined (off-line) by analyzing the charged particle multiplicity measured by two 
stations of silicon strip detectors (MSD) which sample the pseudorapidity intervals 
2 < r] < 3 and 3 < T] < 4. 



Transverse mass spectra of strange particles in Pb-Pb at 158 A GeV/c 



5 



3. Data analysis 

The results presented in this paper are based on the analysis of the full data sample 
consisting of 460 M events of Pb-Pb collisions. 

The strange particle signals are extracted using geometric and kinematic constraints, 
with a method similar to that used in the WA97 experiment As an example of the 
quality of the NA57 data, the invariant mass spectra for 7r + 7r~, p7r, An and AK after 
all analysis cuts are shown in figure El Hyperon and Kg peaks are well centered at the 




0.45 0.5 0.55 0.6 

M(jiV) [GeV] 



1.25 




1.25 



M(pjt") [GeV] 



M(p;i + ) [GeV] 



> : 

CD1200- 
CM 

UU000 - 

CD 
> 
CD 

_ 800- 

o 



600- 




> 300 
CD 



~~ 250 

to 

c 

CD 

a> 200 



150 



100 




1.3 1.35 1.4 1.45 

M(Aji") [GeV] 



> 

CD, 



05 



§140 

> 

CD 



-100 
80 
60 
40 
20 



ST 



1.3 1.35 1.4 1.45 

M(Ajt + ) [GeV] 



1.6 1.7 1.8 1.9 2 

M(AK") [GeV] 




1.61.7 1.8 1.9 

M(AK + ) [GeV] 



Figure 2. Sample invariant mass spectra for it + tt , pn, Air and AK. 

PDG values PH with FWHM of about 10 MeV (15 MeV for the Kg). 

Although the mass spectra show very low background, a detailed study of the 
residual combinatorial background has been performed for the high statistics samples 
(Kg, A and A) using the method of event mixing. Fake A, A or Kg (V^ ) candidates are 
built by pairing all the negative particles from one event with all the positive ones from 



Transverse mass spectra of strange particles in Pb-Pb at 158 A GeV/c 



6 



a different event, selecting events which are close in multiplicity. Then the fake V°s 
from mixed events are processed with reconstruction and analysis programs, like real 
events. In this way we obtain a good description of the combinatorial background. The 
absolute normalization is fixed by the number of pairs of oppositely charged particles 
in real and mixed events. Figure El shows the tt + tt~ invariant mass distribution for real 
and mixed events before (left) and after (right) the application of the analysis cuts. The 




Figure 3. The tt + tt~ invariant mass distribution for real and mixed events before 
analysis cuts (left) and after analysis cuts (right). The inset in the right plot is a zoom 
on the vertical axis. 

total amount of combinatorial background is estimated to be 0.7%, 0.3% and 1.2% for 
Kg, A and A, respectively. We therefore neglect it. Additional details of this method 
can be found in reference jllj . 

The estimates of the S and Q background, evaluated with a similar technique, are less 
than 4% and 6%, respectively. They have also been neglected for this analysis. 

For each particle species we define the fiducial acceptance window using a Monte 
Carlo simulation of the apparatus, in order to exclude the borders where the systematic 
errors are more difficult to evaluate. The acceptance regions are shown in figure EJ For 
all hyperons, the acceptance window for particle and antiparticle are the same, due to 
the reflection symmetry of our apparatus with respect to the magnetic field direction. 
As a further check, the orientation of the magnetic field was periodically inverted. 

All data are corrected for geometrical acceptance and for detector and 
reconstruction inefficiencies on a event by event basis, with the following procedure: 

• for each reconstructed particle a sample of Monte Carlo particles is generated with 
the same measured rapidity y and transverse momentum p T , random azimuthal 
angle and the primary vertex position generated randomly within the target 
according to the measured beam profile; 

• the Monte Carlo particles are traced through a simulation of the apparatus § based 
§ Silicon pixel detector efficiencies are accounted for chip-by-chip. Each chip contains 2032 or 1006 
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Figure 4. The y-pr acceptance windows superimposed to the data samples. Dashed 
lines show the position of mid-rapidity (y cm = 2.92). 



on GEANT [TJ] allowing them to decay according to their life times and random 
internal decay angles. 

• the hits of the decay tracks are merged with those of a real event — with a hit 
multiplicity close to that of the original event — in order to account for background 
tracks and electronic noise; 

• the merged event is processed through the NA57 reconstruction and analysis chain; 

• the weight to be associated with the real event is calculated as the ratio of the 
number of Monte Carlo generated events to the number of Monte Carlo events 
successfully reconstructed and retained by all analysis selection criteria. 

It has been checked that the final weights are not sensitive to the experimental smearing 
in pi and y. 

The simulation used for calculating the correction factors has been checked in detail [T3] 
by comparing real and Monte Carlo distributions for several parameters, and they show 
good agreement. As an example, we show in figure El three such distributions for each 

sensitive cells (pixels) where the pixel size is, respectively, 50 x 500 /im 2 or 75 x 500 /zm 2 [H]. 
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particle species. 

The systematic errors on the inverse slopes have been estimated to be about 10% 
for all the strange particles. 

The weighting method described above is CPU intensive; therefore, while each 
of the reconstructed 5s and fls have been individually weighted, for the much more 
abundant Kg, A and A samples we only weighted a small fraction of the total sample in 
order to reach a statistical accuracy better than the limits imposed by the systematic 
error. Table [T] shows the strange particle samples collected by NA57 and those used for 
the present analysis. 

Table 1. Statistics of the strange particle samples used in this analysis (individually 
weighted) and collected by NA57 in Pb-Pb collisions at 158 A GeV/c. 



Particle 


K° 


A 


A H~ 




Q- 


n + 


weighted 


3340 


2350 


2718 5858 


1522 


432 


192 


collected 


x400 


x400 


x50 xl 


xl 


xl 


xl 



The charged particle multiplicity distribution, which is shown in figure has been 
divided into five centrality classes (0,1,2,3,4), class being the most peripheral and 
class 4 the most central. The drop at low multiplicities is due to the centrality selection 
applied at the trigger level. The fractions of the inelastic cross-section for the five classes, 
calculated assuming a total cross section of 7.26 barn, are given in table El 



Table 2. Centrality ranges for the five classes defined in figure |5| 



Class 
cr/cr inel (%) 



12 3 4 

53-40 40-23 23-11 11-4.5 4.5-0 



The procedure for the measurement of the multiplicity distribution and the 
determination of the collision centrality for each class is described in reference |14j . 

4. Inverse slopes of the transverse mass spectra 

The double-differential (y, m T ) distribution for each particle species has been 
parametrized using the expression 

d 2N x { rn T \ 

assuming the rapidity distribution to be flat within our acceptance region (f(y) = const). 
The inverse slope parameter T app ("apparent temperature") has been extracted by 
means of a maximum likelihood fit of equation 121 to the data. As discussed in the 
next section, this apparent temperature is interpreted as due to the thermal motion 
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Figure 5. Comparison between real (points with error bars) and Monte Carlo (full 
lines) distributions. 1 st row, Kg: the decay length of the Kg 1 ; the closest distance 
in space (close parameter) between the extrapolated tt + and 7r~ tracks coming from 
the Kg decay; the cosine of 9*, the angle between the tt~ and the Kg lines of flight, 
in the Kg reference system. 2 nd row, A+A: the decay length of the A; the close 
parameter between the extrapolated p and ir; the azimuthal decay angle < 



row, 



, — + 



: the decay length of the S; the close of the A coming from the 5; the 7r impact 
parameter. 4 th row, the decay length of the f2; the azimuthal decay angle 

<fi in the Q decay; the azimuthal decay angle <p in the subsequent A decay. 
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Figure 6. Charged particle multiplicity distribution and limits of the five centrality 
classes. For each class the corresponding cross-section is indicated. 



coupled with a collective transverse flow of the fireball components assumed to be in 
thermal equilibrium. 

In the fits of the Q~ and Q particle spectra, the lower limit for the accepted transverse 
momentum was chosen to be 1.4 GeV/c. This choice allows us to exclude a possible 
instrumental bias at low p T for Q. The values of the transverse expansion velocity and 
of the freeze-out temperature obtained within the blast-wave model (see section 5) do 
not depend significantly on the choice of the lower limit. 

The 1 / m T dN/ dm T distributions for strange particles measured for the full centrality 
range under consideration are shown in figure [7| The shapes of all spectra are well 
described by exponential functions. In the next section, we exploit the deviations from 
the exponential in order to disentangle the transverse flow from the thermal motion. The 
inverse slope parameters T app of the transverse mass spectra are given in table El They 
are in agreement within the errors with those measured over a smaller centrality range 
(about the 40% most central inelastic Pb-Pb cross-section) by the WA97 experiment |15j . 



Table 3. Inverse slope parameter T app (MeV) of the strange particles in the full 
centrality range (0-4). The first error is statistical, the second one systematic. 





A 


A 


237 ± 4 ± 24 


289 ± 7 ± 29 


287±6±29 









n + 


297±5±30 


316 ± 11 ±30 


264 ± 19 ± 27 


284 ± 28 ± 27 
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Figure 7. Transverse mass spectra of strange particles for the 53% most central Pb- 
Pb cross-section. The superimposed exponential functions have inverse slopes equal to 
the T app values obtained from the maximum likelihood best fits. In the Cl~ and fl 
spectra, the open points have not been used in the best fit calculation. 



A compilation of data|| on the m T inverse slopes as a function of the particle mass 
in Pb-Pb interactions at 158 A GeV/c is shown in figure |H1 

As already found by the WA97 experiment, the inverse slopes of the Q hyperon and of 
the 5 fall below the line drawn through the n, K and proton points. This observation 
has been interpreted as due to an early freeze-out of multi-strange hadrons [2%] . 

|| NA49 results are taken from the following references: 7r _ , K~ and K + from JSj; Kg from |17|: 7r + , p 
and p from JSj; </> from A and A from |2J]]; S~ and S + from [21]; £l~ and fl + from (221; deuteron 
from [231 - NA44 results are collected for 7r _ , ir + , K~, K + , p, p from reference [21] and for the deuteron 
from reference [2S]- The point at m ~ 3.1 GeV/c 2 corresponds to the NA50 determination of a thermal 
spectrum for the J/i[> charmonium state [20] ■ The J/ip spectra are fitted using a modified Bessel 
function, namely 1/T ■ rn\ ■ K\ (m t /T) ; this yields an inverse slope T = 255 ± 4 MeV (the value showed 
in figure [SJ for the most central collisions, which is smaller by about 5-10 MeV than the one obtained 
with the exponential function [20] • Finally the NA50 inverse slopes for the p + u> and <fi vector-mesons 
are taken from reference |27| . 
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Figure 8. Inverse slopes as a function of the particle rest mass in Pb-Pb interactions 
at 158 A GeV/c (see text for references). Only statistical errors are shown. The line 
through 7r, K and proton is shown to guide the eye. 

4-1. Centrality dependence 

We have fitted the weighted data by equation|2]separately for each of the five multiplicity 
classes indicated in figure E3 In order to reduce the statistical errors, the Q~ and Q + 
samples have been divided into only three centrality classes instead of five: + 1,2 and 
3 + 4. The 1 / m T dN / dm T spectra are shown in figure |H1 together with the exponential 
functions having inverse slopes as obtained from the maximum likelihood fit. The 
numerical values of the inverse slopes are given in table 0] An increase of the inverse 
slopes with the centrality is observed in Pb-Pb for A, H + and possibly also for A. 
Inverse slopes for p-Be and p-Pb collisions as measured by WA97 are also given in 
table El In central and semi-central Pb-Pb collisions (i.e. classes 1 to 4) one observes 
baryon-antibaryon symmetry in the shapes of the spectra: this can be seen in the values 
of the inverse slopes which are compatible for each hyperon and its antiparticle. This 
symmetry is not observed in p-Be collisions. 

The similarity of baryon and antibaryon m T slopes observed in Pb-Pb suggests that 
strange baryons and antibaryons are produced by a similar mechanism in nuclear 
collisions. 
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Figure 9. Transverse mass spectra of the strange particles for the hve centrality classes 
of figure H3 For each species, class 4 is displayed uppermost and class lowermost. In 
the J7 _ +f2 spectra, the open points have not been used in the best fit calculation. 



Table 4. Inverse m-i slopes (MeV) of the strange particle distributions for the five Pb- 
Pb centrality classes, and for p-Be and p-Pb interactions [2H1- Only statistical errors 
are shown. In Pb-Pb, systematic errors are assumed to be independent of centrality. 





Pb-Pb 




p-Be 


p-Pb 





1 


2 


3 


4 




197 ±4 


217 ±6 


239 ± 15 


239 ±8 


233 ±7 


244 ±8 


234 ±9 


A 


180 ±2 


196 ±6 


237 ± 19 


274 ± 13 


282 ± 12 


315 ± 14 


305 ± 15 


A 


157 ±2 


183 ±11 


277 ± 19 


264 ± 11 


283 ± 10 


313 ± 14 


295 ± 14 




202 ± 13 


235 ± 14 


290 ± 20 


290 ± 11 


295 ±9 


304 ± 11 


299 ± 12 




182 ±17 


224 ± 21 


232 ± 29 


311 ±23 


294 ± 18 


346 ± 28 


356 ± 31 




169 ± 40 


334 ± 99 


274 ± 34 


274 ± 28 


268 


±23 



5. Blast-wave description of the spectra 



In this section we use the statistical hadronization model of reference J5 to describe the 
strange particle spectra discussed above. The model assumes that particles decouple 
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from a system in local thermal equilibrium with temperature T, which expands both 
longitudinally and in the transverse direction; the longitudinal expansion is taken to 
be boost-invariant and the transverse expansion is defined in terms of a transverse 
velocity field. Finally, the statistical distributions are approximated by the Boltzmann 
distribution. 

The blast-wave model predicts a double differential cross-section for a particle species j 
of the form: 

d 2 N~ , f R ° /m T coshp\ fp T smhp\ , . . 

where p(r) = tanh -1 /3±(r), K\ and 7 are two modified Bessel functions and Aj is a 
normalization constant. With respect to a cylindrical reference system (r,<p,z,t), the 
freeze-out hypersurface is constrained by < r < Rq, < <p < 2tt and dtf/dr = 0; the 
last condition assumes that the particles decouple suddenly from the whole transverse 
profile of the fireball at time tf. In these expressions Rg is the transverse geometric 
radius of the source. 

In case of a peripheral collision, the azimuthal symmetry is evidently broken; however 
it is recovered when considering the m T spectrum of particles accumulated over many 
events with random impact parameters, which is our approach. 
The transverse velocity field /3±(r) can be parametrized according to a power law: 



r < R G (4) 



Rg. 

With this type of profile the numerical value of Rg does not influence the shape of the 
spectra but just the absolute normalization (i.e. the Aj constant). Once the transverse 
flow profile (i.e. equationHJ) has been fixed, the shape of each spectrum is determined by 
the temperature, the velocity of the transverse expansion and the mass of the particle. 
The parameters which can be extracted from a fit of equation |3] to the experimental 
spectra are thus the thermal freeze-out temperature T and the surface transverse flow 
velocity /3s- In order to compare results from different profile hypotheses, corresponding 
to different values of the exponent n in equation 0J the average transverse flow velocity 
can be used instead. Assuming a uniform particle density, the average^ is |3Hj : 

5.1. Global fits to the spectra using different flow velocity profiles 

Preliminary NA57 results of the blast-wave model were presented using a constant 
(n = 0) velocity profile, for the purpose of comparison with the NA49 analysis [32] . In 

% A more sophisticated averaging can be achieved by incorporating not only the transverse geometry 
of the model but also the phase space density of particles [3D]. According to this definition, < j3±_ > is 
also a function of T and it differs from the values calculated according to equation[S]by 2% if n = 1/2, 
by 5% if n = 1 and by 10% if n = 2; obviously for n — 0, < /3j_ >= Ps independently of the average 
definition. 
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the following, we will adopt primarily the linear (n = 1) r-dependence of the transverse 
flow velocity. In fact, a constant flow profile description of the transverse expansion is 
internally inconsistent, because the outermost fluid elements of the expanding fireball 
should have travelled faster than the innermost ones: (3{r\) > /3(r 2 ) if r\ > r 2 . On the 
other hand, the quadratic profile (n = 2) is disfavoured by data, as will be shown below. 
We also consider the n = 1/2 profile which has been suggested to resemble closely the 
solution of the full hydro-dynamical calculation |33| . 

The global fit of equation|S]with n = 1 to the data points of all the measured strange 
particle spectra is shown in figure [TUl it successfully describes all the distributions with 
X 2 / n df = 37.2/48, yielding the following values for the two parameters T and < (3± >: 

T = 144 ± 7(stat) ± 14(syst)MeV , < (3±>= 0.381 ± 0.013(stat) ± 0.012(syst) . 




m T - m [GeV] <P ± > 

Figure 10. Blast- wave fit to the spectra with a linear profile (left); corresponding la 
and 3cr confidence level contours in the < (3± >-T plane (right). 

The T and < (3± > parameters are found to be statistically anti-correlated, as can be 
seen from the confidence level contours shown in figure Conversely, the systematic 
errors on T and < (3± > are found to be correlated and they are estimated to be 10% 
and 3%, respectively. 

The results of the fits with all the four profile hypotheses are given in table 03 
The three profiles n = 0, n = 1/2 and n = 1 yield similar values of the freeze-out 
temperatures and of the average transverse flow velocities, with good values of x 2 /ndf . 
The quadratic profile is disfavoured by our data. 

For m T sufficiently larger than the particle mass mo, the inverse slope T app is 
expected to be independent of m jOJ- In the non-relativistic regime (p T < m ) T app 
is expected to increase with the rest mass of the particle; e.g. in reference [SI] it has 
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Table 5. Results of the blast-wave model fit using different transverse flow profiles. 
The quoted errors are statistical. The systematic errors on the temperature and on the 
velocities are estimated to be about 10% and 3%, respectively, for all the four profiles. 





n = 


n = 1/2 


n — 1 


n = 2 


T (MeV) 


158 ±6 


152 ±6 


144 ± 7 


151 ± 11 


Pa 


0.396 ±0.015 


0.493 ±0.016 


0.571 ±0.019 


0.633 ±0.028 


<P±> 


0.396 ±0.015 


0.394 ±0.013 


0.381 ±0.013 


0.316 ±0.014 


X 2 /ndf 


39.6/48 


36.9/48 


37.2/48 


68.0/48 



been derived that T app « T + |m < /3± > 2 : the larger the mass of the particle the more 
pronounced the flattening of the spectra at low m T . The deviations from exponential 
behaviour are larger for large flows and at low m T . Due to its light mass, the Kg 
spectrum is predicted by the model, and indeed observed, to be steeper at low m T . Such 
behaviour is illustrated in figure ITT1 where we compute the "local" inverse slope, i.e. the 
inverse slope T app of the exponential function exp(— m T /T app ) which is tangent to the 
blast-wave curve, as a function of m T — mo- Pions, which in principle would be an even 
better probe for this behaviour, are unfortunately affected by resonance decays (at low 
p T ) which cannot be neglected |3SHSE!- 




Figure 11. Local inverse slopes according to the blast-wave model using the n = 1 
transverse profile. 

The particles have been divided into two groups — those which share valence quarks 
with the nucleons and those which do not — since it is known that the particles of 
the two groups may exhibit different production features. The fit procedure has been 
repeated separately for the two groups (n = 1 case) and the results are given in table El 
They suggest common freeze-out conditions for the two groups. Since the interaction 
cross-sections for the particles of the two groups are quite different, this finding would 
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Table 6. Thermal freeze-out temperature (T), surface {(is) and average (< j3±_ >) 
transverse flow velocity in the full centrality range, assuming a linear transverse profile 
(n = f ). The first error is statistical, the second one systematic. 



particles 


T (MeV) 


Ps 


<P±> 


X 2 /ndf 


K° s , A, S- 


146 ± 8 ± 14 


0.564 ±0.023 ±0.018 


0.376 ±0.015 ±0.012 


18.1/23 


A, H + , n~, Q + 


130 ± 28 ± 14 


0.604 ±0.048 ±0.018 


0.403 ±0.032 ±0.012 


18.5/23 



suggest limited importance of final state interactions (i.e. a rapid thermal freeze-out) 
and a similar production mechanism. A similar conclusion concerning the evolution of 
the system was reached from the study of the HBT correlation functions of negative 

pions ma main- 

It has been suggested [2H], based on WA97 results on the hyperon m T slopes 
compatible with those discussed in section 4, that the thermal freeze-out occurs earlier 
for Q and possibly for H than for particles of strangeness 0, due to the low scattering 
cross-sections for Q and S. The la contours of the separate blast wave fits for singly and 
multiply strange particles are shown in figure^! Both groups of particles are compatible 
with the global-fit determination. However, the fit for the multiply strange particles is 
statistically dominated by the S; in fact the 5±f2 contour remains essentially unchanged 
when fitting the H alone. Therefore we can only conclude that the H undergoes a thermal 
freeze-out which is compatible with that of A and Kg. 

Due to the lower statistics of the Q, from its spectrum alone it was not possible to 
extract significant values for the two freeze-out parameters, as can be done for the 5. 
Any possible deviation of the Q from the observed common freeze-out undergone by Kg, 
A and H can only be inferred from the integrated information of the Q spectrum, i.e. 
from its inverse slope. 

In figure HH1 we plot the inverse slopes of figure |H1 superimposing the blast-wave 
model results. The full lines represent the blast-wave calculations when fitting the m T 
spectrum of a particle, as a function of mass, in the common range 0.05 < m T —m < 1.50 
GeV/c 2 , for two different freeze-out conditions: absence of transverse flow (< (3± >= 0) 
and our best fit determination. The inverse slope is also a function of the m T — mo 
range where the fit is performed. This dependence is stronger for heavier particles (see 
figure ITTjl . Therefore we have also computed the blast- wave inverse slopes of Kg, A, H 
and Q spectra in the m T — m ranges of NA57 (closed circles in figure H*3*j) . The measured 
inverse slopes of Q show a deviation from the global trend of the other strange particles. 

5.2. Centrality dependence 

We have performed the global fit to the spectra for each of the five centrality classes 
defined in section 3. In figure El we show the la confidence level contours in the 
< p± >-T plane as obtained for the n = 1 profile. The numerical values of the fit 
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!—< 160 
> 




Figure 12. The thermal freeze-out temperature vs average transverse flow for the 
blast-wave fits (n — 1 profile) to multiply and singly strange particle spectra. The la 
contours are shown, with the marker indicating the optimal fit location. The global fit 
contour (all) is also shown. 




Figure 13. Prediction of the blast-wave model for inverse slopes (see text for details). 
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parameters are given in table [H The observed trend is as follows: the more central 
the collisions the larger the transverse collective flow and the lower the final thermal 
freeze-out temperature. The higher freeze-out temperature for peripheral collisions may 
be interpreted as the result of an earlier decoupling of the expanding system. 
A similar centrality dependence of the freeze-out parameters has been observed by the 
PHENIX experiment in Au-Au collisions at ^Jsnn — 130 GeV from the blast-wave 
analysis of 7r ± , K ± , p and p spectra |4~U] . In particular, PHENIX has measured! 
T = 121 ± 4 MeV, < /3j_ >= 0.47 ± 0.01 with x 2 jndf = 34.0/40 for the most 
central 5% of the Au-Au cross-section and T = 140 ± 4 MeV, < (3 X >= 0.39 ± 0.01 
with x 2 / n df — 68.9/40 for peripheral collisions (30-60% centrality, which corresponds 
roughly to the NA57 class 0). For the same colliding system, the STAR experiment has 
also determined |H] the thermal freeze-out for central collisions (see table |7|): STAR 
results would suggest a rather earlier decoupling of the multi-strange H with respect to 
non strange and singly strange particles at RHIC energy. 

We also quote the result of our fit for the most central class, when the n = profile is 
assumed: T = 134±10±13 MeV, < p ± >= 0.457±0.025±0.014 and x 2 /ndf = 47.2/43. 
This result is found to be consistent within two standard deviations with that from a 
similar analysis performed by NA49§ assuming the constant profile hypothesis [32J. 
In table [7| we also make a comparison with results from different measurements which 
exploit the p r dependence of the HBT correlation functions of negative pions |SZ1 EHl EHj • 
The thermal freeze-out temperatures and transverse flow velocities are found to be in 

% In the PHENIX analysis, statistical and systematic errors are added in quadrature before the fit. 
§ The NA49 centrality ranges are 0-5% for K± and (f>, 0-10% for tt, A and S, 0-20% for fi. 
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Table 7. Thermal freeze-out temperature and transverse flow velocity in the five 
centrality classes, assuming a linear transverse profile (n = 1). Only statistical errors 
are quoted; the systematic errors on T and < (3± > are estimated to be 10% and 3%, 
respectively. We also give the blast- wave determinations of T and < /3± > by the NA49 
experiment in Pb-Pb collisions at 158 A GeV/c [32], and by the PHENIX and STAR 
experiments in Au-Au collisions at ^/snn = 130 GeV ^Ul^D- Finally, we quote the 
measurements by the NA49 [37], WA97 [38] and CERES [3J| experiments which use 
HBT correlation combined with information from h~ and/or deuterium single particle 
spectra. 



Exp. 






technique 


T (MeV) 


< 0± > 


X 2 /ndf 


NA57 


40-53 


17.3 


this analysis, n = 1 


139 ± 25 


0.35 ±0.04 


33.2/35 


NA57 


23-40 


17.3 


this analysis, n = 1 


151 ± 15 


0.35 ±0.03 


42.3/42 


NA57 


11-23 


17.3 


this analysis, n = 1 


146 ± 17 


0.37 ± 0.03 


41.3/43 


NA57 


4.5-11 


17.3 


this analysis, n = 1 


135 ± 14 


0.41 ±0.02 


31.8/41 


NA57 


0-4.5 


17.3 


this analysis, n — 1 


118 ± 13 


0.45 ± 0.02 


52.7/43 


NA49 


0-5,10,20 


17.3 


K+,p,A,£-,Q-, n = 


127 ±4 


0.48 ± 0.01 


118/43 


NA49 


0-5,10,20 


17.3 


K-,p,M,S + ,n + , n = 


122 ±2 


0.48 ± 0.01 


46/41 


PHENIX 


60-92 


130 


7i' ± ,K ± ,p,p, n = 1 


161±g 




36.3/40 


PHENIX 


30-60 


130 


7i" ± ,K ± ,p,p, n = 1 


140 ±4 


0.39 ±0.01 


68.9/40 


PHENIX 


15-30 


130 


7i" ± ,K ± ,p,p, n = 1 


134 ±2 


0.43 ±0.01 


36.2/40 


PHENIX 


5-15 


130 


7i" ± ,K ± ,p,p, n = 1 


125 ±2 


0.46 ±0.01 


34.7/40 


PHENIX 


0-5 


130 


7i" ± ,K ± ,p,p, n = 1 


121 ±4 


0.47 ±0.01 


34.0/40 


STAR 


0-10 


130 


S- S + , n= 1/2 


182 ±29 


0.42 ±0.06 


13/15 


STAR 


0-10 


130 


7r,K,p,A, n = 1/2 


103 ±7 


0.57 ±0.01 




NA49 


0-3 


17.3 


HBT + h~ + d 


120 ±12 


0.55 ±0.12 




WA97 


25-40 


17.3 


HBT + br 


140+g 


30+ 09 




WA97 


12-25 


17.3 


HBT + br 


121+15 


47+ - 07 
u -^' -0.10 




WA97 


5-12 


17.3 


HBT + Ti- 


n7±i? 


0.48lg-$S 




WA97 


0-5 


17.3 


bet + h~ 


120115 


46+ 07 
u -* D -o.io 




CERES 


15-30 


17.3 


HBT 


fixed at 120 


u -* y -0.06 




CERES 


10-15 


17.3 


HBT 


fixed at 120 


u -°°-0.05 




CERES 


5-10 


17.3 


HBT 


fixed at 120 


46+ 04 




CERES 


0-5 


17.3 


HBT 


fixed at 120 


u - oo -0.03 





good overall agreement with those extracted from the blast-wave fits. 
6. Conclusions 

We have analysed the transverse mass spectra of high statistics, high purity samples 
of Kg, A, H and fl particles produced in Pb-Pb collisions at 158 A GeV/c over a wide 
range of collision centrality (the most central 53% of the Pb-Pb inelastic cross-section). 

The inverse slopes agree with those measured by WA97. For each hyperon species, 
the m T slope is found to be in good agreement with that of the antiparticle. An increase 
of the inverse slopes with the collision centrality is observed for A and S + . 

The analysis of the transverse mass spectra at 158 A GeV/c in the framework of the 
blast-wave model suggests that after a central collision the system expands explosively; 
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the system then freezes-out when the temperature is of the order of 120 MeV with 
an average transverse flow velocity of about one half of the speed of light. Particles 
with and without valence quarks in common with the nucleon appear to have a similar 
behaviour. The inverse slope of the Q particle deviates from the prediction of the blast- 
wave model tuned on the other strange particles (Kg, A and 5). Thermal freeze-out 
conditions (i.e. temperature and transverse flow velocity) are found to be similar to those 
measured at RHIC energy for singly strange and non strange particles, in both central 
and peripheral collisions; conversely the temperature measured for the multiply strange 
S is significantly larger at RHIC. Finally, the results on the centrality dependence of the 
expansion dynamics indicate that with increasing centrality the transverse flow velocity 
increases and the final temperature decreases. 
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